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Classically it is impossible to have transport without transit, i.e., if the points one, two and 
three lie sequentially along a path then an object moving from one to three must, at some point in 
time, be located at two. However, for a quantum particle in a three- well system it is possible to 
transport the particle between wells one and three such that the probability of finding it at any time 
in the classically accessible state in well two is negligible. We consider theoretically the analogous 
scenario for a Bose-Einstein condensate confined within a three well system. In particular, we 
predict the adiabatic transportation of an interacting Bose-Einstein condensate of 2000 ^Li atoms 
from well one to well three without transiting the allowed intermediate region. To an observer of this 
macroscopic quantum effect it would appear that, over a timescale of the order of Is, the condensate 
had transported, but not transited, a macroscopic distance of ^ 20/im between wells one and three. 



The system under consideration is schematically shown 
in Fig. 1(a), where a three-dimensional harmonic trap is 
split into three regions via the addition of two parallel re- 
pulsive Gaussian potentials. With the Bose-Einstein con- 
densate (BEG) [blue object in Fig. 1(a)], initially in well 
1, we show how it is possible, through adiabatic changes 
to the tunneling rates between the wells, to transport 
it into well 3 with minimal (ideally zero) occupation of 
the intervening well. This effect as a function of time is 
shown in Fig. 1(b), where an interacting BEG of 2000 
^Li atoms is transported from well 1 to well 3 over a 
timescale of ~ Is, with less than 1% atoms occupying 
well 2 at any particular time. As such it appears that 
the BEG is transported from well 1 to well 3 without 
transiting through well 2. 

This effect of transport without transit (TWT) can 
be likened to the lay concept of teleportation. However, 
although TWT relies on quantum control of the global 
BEG state and associated tunneling matrix elements, it 
is quite distinct from the quantum definition of telepor- 
tation [1 . In the TWT of a BEG we describe the many 
body system in a time dependent mean-field approxima- 
tion. As such the wavefunction used to describe the con- 
densed state is a classical field and can not describe such 
properties as entanglement and hence quantum telepor- 
tation. 

The ideas underpinning the protocol for TWT stem 
from Stimulated Raman Adiabatic Passage (STIRAP) 
[2j |3l m [5]. STIRAP is a robust optical technique for 
transferring population between two atomic states, |1) 
and 1 3), via an intermediate excited state, |2). Using 
off-resonant pulses to couple states |1) to |2) and |2) to 
1 3), characterised by coupling parameters K12 and i^23, 
and such that K23 precedes and overlaps i^i2, the pop- 
ulation can be adiabatically transferred from state |1) to 
1 3). Population transfer is achieved via a superposition of 
states |1) and |3) with the occupation of state |2) strongly 
suppressed. These techniques are used in quantum op- 
tics for coherent internal state transfer [H El |7l [8] and 
have been proposed for applications in three channel opti- 
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FIG. 1: Transport without transit (TWT) of a BEC in a three 
well system, (a) Schematic representation of our system at 
t=0 in the {z — y) plane. Two parallel, repulsive Gaussian 
barriers embedded in an ambient harmonic trap divide the 
system into three wells, with the BEC initially occupying well 
one. (b) Isosurface plots of atomic density (riiso = O.lno, 
where no is the initial peak density of the BEC) showing the 
adiabatic transportation of a BEC of 2000 ^Li atoms over a 
distance of 20/im (simulated using the 3D Gross-Pitaevskii 
equation) . 

cal waveguides [9^ while an analogous approach has been 
proposed for state transfer from one atom laser beam to 
another [10]. Recently this protocol has been proposed 
to transport single atoms [HI HI], Gooper pairs [13 and 
electrons [14, 15 , 16 , 17^. Here we extend these ideas to 
the transport of dilute gas BEGs containing thousands of 
atoms. 

In this article we elucidate the properties of the three- 
well system by first considering a three-mode approxi- 
mation [iSl [191 EQ], where the form of the potential is 



2 



not important. We then employ the mean-field Gross- 
Pit aevskii equation (GPE) to qualitatively describe the 
BEG dynamics and consider experimental scenarios in 
which to realise macroscopic matter- wave TWT. 

Reducing our three- well system, shown in Fig. 1(a), 
such that each well is described by a single mode basis 
[IHl [m [20] , ^i, enables its properties to be described via 

^{t) = Mt)^i + Mt)^2 + Mt)'^s (1) 



where 





Ui -Ki2 
-Ki2 U2 —K23 
-K23 Us 



The amplitude of each mode is expressed as i/ji = 
y/Nie^^\ where Ni and 6i are the occupation and phase 
of the ith mode, respectively {i = 1,2,3). The system 
is normalised such that Xli=i^i(^) — where Nt 
is the total number of atoms in the system. The pa- 
rameters K12 and K2S describe the wavefunction over- 
lap, and hence tunneling rate, between wells 1 and 2, 
and 2 and 3, respectively. Furthermore, the dimension- 
less on-site interaction energy per particle is given by 
Ui = ^QmNi/NT^ where is the groundstate energy 
of well i and Qm is a dimensionless parameter describing 
the nonlinear atomic interactions within the system . 

The modulation of the wavefunction overlaps K12 and 
K2S controls a transfer of atoms between the wells. We 
assume these parameters vary with time as. 



Ki2{t) 



[7rt/(2tp)], 
■[7rt/(2tp)] 



(3) 



where tp is the total pulse time and the maximum tunnel- 
ing rate is defined by l^max- We employ this protocol due 
to it robustness against non-linear effects arising from the 
inter- atomic interactions at t = and t = tn. 



In the limit tp ^ 00 and for gr, 
the modes are given by ^ , 



the evolution of 



= sin 61 sin 62 ^1 + cos 62 ^2 + cos 61 sin 62 ^3 
D_ — sin Oi cos O2 ^1 — sin 02 ^2 + cos 01 cos 02 ^3 
Do cos 01 ^1 - sin 0i ^3 (4) 



where 



01 = arctan [K12/K23] 
1 
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The corresponding mode-energies are shown in Fig. 2(a). 
For an initial state where all of the atoms are in well 1 
the system adiabatically follows the green/middle line. 
This corresponds to the passage of atoms from well 1 to 
well 3, with a heavily suppressed occupation of well 2, as 
shown in Fig. 2(b), corresponding to TWT. 
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FIG. 2: Dynamics of the system according to the three-mode 
analysis, (a) Energies of the eigenstates Do and D- 

of the non-interacting (gm = 0) system, (b) Evolution of 
Ni{t)/NT and N3{t)/NT for g^ = and tp oo. (c) 
Ns{t = tp)/NT as a function of gm and tp, with white and 
black representing N3(t — tp)/NT = 1 and respectively, (d) 
max[A^2(t)/A^T] as a function of gm and tp, with white and 
black representing A^3(t = tp)/NT = and 1 respectively. In 
(c) and (d) the solid black and blue/grey curves represent 
Nsit = tp)/NT = 0.99 and max[N2{t)/NT] = 0.01 respec- 
tively, and the region bounded by both corresponds to high 
fidelities e > 0.99. We have assumed = = and 
E2 = O-l/i^max throughout. 

Figures 2(a,b) are in an ideal limit where the atomic 
interactions are zero {gm = 0) and the time over which 
the pulses were applied was large {tp oo). However, 
for a realistic system it is instructive to examine how this 
ideal picture changes as the pulse times (tp) are reduced 
and interactions are included [18^. It is possible to pa- 
rameterise the evolution of the system via the fraction of 
atoms in well 3 at the end of the protocol, Ns{tp)/NT, 
and the maximum number of atoms occupying well 2 
during the protocol, mdix[N2{t) /Nt]. These quantities 
are mapped out in Figs. 2(c) and (d), respectively, as a 
function of the strength of the nonlinear interactions gm 
and the pulse time tp, with efficient TWT occurring in 
the white regions. Defining a fidelity, e, for TWT through 
Ns{tp)/NT > e and max[7V2(^)/A^T] < 1 - e we find that 
to achieve e > 0.99 we require that l^^nl < ^2 ~ 
^p^max^400. In the absence of nonlinear interactions 
(^ = 0) the condition for tpl^max > 400 comes from the 
adiabatic limit of the system and is governed by the en- 
ergy difference between the groundstate energies of the 
wells. As noted by Graefe et al. [18 , the introduction of 
nonlinear interactions introduces new nonlinear "eigen- 
states", which can inhibit adiabatic transfer. We note 
that in their approach they considered a Gaussian tun- 
neling scheme. However, the protocol which we employ 
[Eq. ([3|] is much more robust to non-linear effects, since 
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the energies of the additional non-hnear states are not 
close to the dark state mediating the transfer. 

The mode analysis presented above gives a qualitative 
description of adiabatic transport for a three well system. 
To investigate TWT quantitatively for realistic scenarios 
the GPE is employed. The GPE mean-field model has 
had great success in describing the dynamics of BECs, 
e.g. the formation of vortex lattices in rotating BECs [22l 
[23] , the quantum reflection of BECs off silicon surfaces 
[24] |25] , the breakdown of Bloch oscillations of BECs 
in optical lattices [26] |27] and the formation of bright 
solitary waves in attractive BECs [2B • 

The GPE model is valid in the limit of zero tem- 
perature and describes the BEG by a macroscopic or- 
der parameter, or "wavefunction" , ^(r, t) which repre- 
sents the mean-field of the Bose-condensed atoms. This 
macroscopic wavefunction can be expressed as ^(r, t) = 
^/n(J~t) ex.-p[iO{r^t)]^ where n(r, t) is the atomic density 
and 6>(r, t) is a macroscopic phase. The evolution of the 
wavefunction ^(r, t) is described by the GPE, 



10' 



ih 



g^(r,t) 
dt 



^(r,t). 



(7) 



Here the nonlinear coefficient is given by ^ = 47r fi^a/m, 
where a is the s-wave scattering length that characterises 
the atomic interactions in the BEG. We assume a trap- 
ping potential F(r, t) of the form. 



m 



Vi2{t) exp 



\2" 



2^21 



2a2 



V23{t) exp 



2a2 



(8) 



The first term defines the cylindrically-symmetric 
parabolic trap with radial and axial trap frequencies 
uj± and ujz^ respectively. The second and third terms 
represent the repulsive Gaussian barriers, positioned at 
z = ±2^0, with width a and time-dependent amplitudes 
Vi2{t) and V23{t). Such barrier potentials can be in- 
duced by the optical dipole force from two parallel blue- 
detuned laser beams [29] or using magnetic fields on an 
atom chip [30]. The barrier amplitudes, which can be 
varied by modifying the laser intensity or the atom chip 
currents, controls the tunnelling rates between the neigh- 
boring wells. 

Recall that in the mode analysis we employed tun- 
nelling rates which initially had opposing values [Eq. (|3|]. 
Due to difficulties in initialising such a state in an ex- 
periment and simulation, we begin our simulations with 
barriers of identical height Knax, and therefore initially 
identical tunnelling rates. The exact functional form of 
Vi2{t) and V23{t) that we employ can be found in the 
Methods section. Over time, the barriers are smoothly 
lowered to a minimum value Vmin before being increased 
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FIG. 3: Dynamics of the system according to the ID GPE. 
(a) Carpet plot showing the evolution of condensate den- 
sity (dark=high density, light=low density) for an effective 
ID interaction parameter of qid — 0.31 and a pulse time of 
tp = lOOOtho- (b) Same as (a) but for a reduced pulse time 
of tp = 14tho- Note the breakdown of the adiabatic trans- 
fer, (c) N3{t = tp)/NT as a function of giu and tp, with 
white and black representing A^3(t = tp)/NT = 1 and re- 
spectively, (d) max[iV2(t)/A^T] as a function of giu and tp, 
with white and black representing ma,'x[N2(t) /Nt] = and 1 
respectively. In (c) and (d) the solid black and blue curves 
represent Nsit = tp)/NT = 0.95 and ma^[N2{t)/NT] = 0.05 
respectively, and the region bounded by both corresponds to 
high fidelities e > 0.95. (e) The evolution of the first nine 
energy eigenvalues (en) as a function of time (^id = 0), with 
the blue, green and red curves corresponding to the D-, Dq 
and D+ mode energies respectively. In all cases Vmin = ^hcoz , 
Vmax = lO^huJz, cr — 0.16/ho and zo = 0.48/ho- 

back to Fmax, such that the tunnelling rate variation ap- 
proximates a Gaussian. Importantly, there is a time de- 
lay of size r between the pulsing of Vi2{t) and V2s{t). 
When Vi2{t) is pulsed before V23(t), we term this the in- 
tuitive protocol, and when Vi2{t) is pulsed after F23(t), 
we term this the counter- intuitive protocol. Only the 
latter protocol is capable of producing TWT and so we 
concentrate on this. 

The difference between the pulsing schemes used for 
the mode analysis and GPE simulations is seen by com- 
paring the adiabatic evolution of the eigenstates in the 
two approaches. Figs. 2(a) and 3(e). Despite the qualita- 
tively different TWT protocols for the mode analysis and 
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GPE simulations, we see qualitative agreement in the re- 
gions of high fidelity. This suggests that the method of 
transporting BECs is not particularly dependent on the 
exact form of Vi2{t) [Ki2{t)] and V23{t) [i^23(^)], as ex- 
pected for an adiabatic protocol. This has been verified 
through the study of several different functional forms 
for Vi2{t) [Ki2{t)] and ¥23(1) [i^23(^)] which ah produce 
qualitatively similar results. 

The one-dimensional equivalent of the GPE can be 
solved numerically with relative ease and so we consider 
this limit first. Employing harmonic oscillator units, as 
outlined in the Methods section, we consider a system 
defined by a = O.lG/ho, ^0 = 3<j and r = tp/10. 

Figures 3(a,b) show the evolution of the condensate 
density 1^(2:) p for different time pulses but the same 
remaining parameters (^id = 0.31, Vmin = ^huj^ and 
^max = lO'^huJz)' In Fig. 3(a) a large time pulse of tp = 
lOOOtho leads to efficient TWT, with the BEG moving 
smoothly from well 1 to well 3 with a minimal occupation 
of well 2. In Fig. 3(b), however, a significantly reduced 
pulse time of tp = 14tho breaks the adiabaticity of the 
process and causes inefficient transfer, with a significant 
population in well two. 

In Figs. 3(c,d) Ns{tp)/NT and mdix[N2{t)/NT] are 
plotted as a function of tp and ^id, in analogy to the 
three- mode results in Figs. 3(c,d). The curves in 
Figs. 3(c) and (d) denote Ns{t'P)/NT = e (black) and 
mdi-x[N2{t) / Nt] = 1 — e (blue) for e = 0.95, according to 
the GPE. These results are qualitatively similar to the 
mode analysis and show a large region of the parameter 
space where efficient TWT can occur. 

We now consider the possibility of producing efficient 
TWT in a realistic BEG system. We performed simula- 
tions of the full 3D GPE. Since strong nonlinear inter- 
actions suppress TWT we focus on a system with weak 
interaction strength, i.e. a small s-wave scattering length 
a and low atom number Nt. Our simulations are based 
on recent ^Li soliton experiments [3T1[32]. These exper- 
iments have two key advantageous features. Firstly, the 
experiments worked with low atom number, with typ- 
ically several thousand atoms in the condensate. Sec- 
ondly, the experiments employed a Feshbach resonance 
to control the s-wave scattering length and, indeed, this 
allowed the use of a low attractive scattering length of 
the order of a = — O.lnm. This means that in principle 
similar experiments could probe high fidelity parameter 
space of g and tp. 

We consider Nt = 2000 and a = — 0.2nm, and real- 
istic parameters for our trapping system: 00^ = = 
27T X 40Hz, a = 1/im, Fmax = lOOhco^ Vmin = Sfio; and 
zq = 3/im. Initially we consider the condensate dynamics 
under the counter-intuitive protocol and for a pulse time 
of tp = 400a;~^ = 1.6s and a pulse delay of r = 0.16s. 
These dynamics are presented in Fig. 1(b) as chrono- 
logical frames of an isosurface of the BEG density and 
in Fig. 5(a) as a carpet plot of the radially-integrated 
axial density. For these realistic parameters we clearly 
see efficient TWT, i.e. the 2000 ^Li atoms are adiabati- 
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FIG. 4: Demonstration of TWT for a realistic EEC of 2000 
^Li atoms and attractive interactions a = — 0.2nm. (a) Car- 
pet plot showing the evolution of the radially-integrated axial 
density for the counter-intuitive protocol (defined by Eqs. ([9| 
and llO}). (b) Same as (a) but for the intuitive protocol 
[V23(tj-^ V^i2(t) and V23{t) Vi2{t)] giving an appreciable 
population in the middle middle well. In each of the plots 
the horizontal dashed lines correspond to the center of the 
Gaussian barriers ±2:0. We assume cj = 27r x 40Hz, a = l/xm, 
zo — 3/im, Vmax — lOOhoo and Vmin — ^hoo. 

cally transported a distance of approximately 20/im with 
negligible occupation of well 2. Grucially, the timescale 
for this process is just under 2s, which is the lifetime 
of such condensates [32]. These results have a fidelity 
e = 0.985, which is limited by the maximum occupation 
of well 2 during the transfer. Up until now we have de- 
fined well 2 as the spatial region [— ^o,^o] for simplicity. 
However, it may be more appropriate to define well 2 
as the classically- allowed region, i.e. the region between 
—zq and zq where the chemical potential of the initial 
state is less than y(r, t). Under this definition we find 
that the maximum atom number at any given time that 
occupies the classically- allowed region is less than 1% of 
Nt^ giving a fidelity of e > 0.99. 

We have also simulated the dynamics of this system 
for the intuitive protocol. Recall that this corresponds 
to when the first barrier is pulsed before the second bar- 
rier. These dynamics are presented in Fig. 5(b), which 
shows the evolution of the radially-integrated axial den- 
sity. Under this protocol we clearly see the macroscopic 
occupation of well 2 during the transfer. Indeed, at a 
single time during these dynamics over 15% of the atoms 
reside in the classically allowed region in well 2, two or- 
ders of magnitude larger than for the counter-intuitive 
protocol. This demonstrates that a straightforward ex- 
perimental confirmation of TWT is to reverse the pulses 
and compare the condensate density in the middle well 
half way through the pulse sequence. 
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In conclusion we propose a novel protocol for the trans- 
port of BECs in three- well systems. This protocol enables 
the adiabatic transport a macroscopic BEC such that the 
transient occupation of the intermediate well is heavily 
suppressed: transport without transit. In particular, we 
have shown that this works within both a three-mode ap- 
proach and a meanfield approximation, where all of the 
modes of the mean-field system are considered. We have 
mapped out the parameter space for which we expect 
transport without transit to occur. Specifically, we have 
demonstrated the transport-without-transit of an inter- 
acting BEC of 2000 ^Li atoms a macroscopic distance of 
20/im over a timescale of 1.8s. This phenomenon is not 
only of interest from the view point of testing the wave 
nature of a dilute gas Bose-Einstein condensate, but also 
paves the way for a new method of control in atom op- 
tical devices. Future extensions to this work include the 
examination of non-meanfield effects, such as quantum 
fluctuations [33 , and the consideration for systems with 
more than three wells [34] . 

The authors acknowledge useful discussions with Si- 
mon Devitt and David Jamieson. This work is funded 
by the Australian Research Council. Additionally JHC, 
ADG and LCLH are supported by the US National Se- 
curity Agency (NSA), Advanced Research and Develop- 
ment Activity (ARDA) and the Army Research Office 
(ARO) under Contract Nos. W911NF-04-1-0290. 



I. METHODS 

To induce an approximately Gaussian modulation of 
the tunnelling rates [15, 18] we need a functional form 
for the barrier heights of. 



In order to ensure that the transport across the three- 
well system is dominated by tunneling, ]4iin is greater 
than the chemical potential of the BEC. 

The dynamics of the BEC are determined through nu- 
merical solutions of Eq. ^ by employing the Crank- 
Nicholson method [35] to evolve the equation. The initial 
state at t = is obtained via propagation in imaginary 
time [36 subject to the number of atoms and the poten- 
tial of Eq. ([8| . The true ground state for the system at 
t = consists of a symmetric state with half the atoms 
in well 1 and half in well 3. To initialize the system we 
set ^(r) =0 for z > 0, thus leaving only atoms in well 
1, as schematically shown in Fig. 1(a). Experimentally 
preferential loading of well 1 can be obtained by an initial 
shift in the parabolic potential [37] . 

For our analysis using the ID GPE we have recast Eq. 
^ in dimensionless form, in terms of harmonic oscillator 
units. The fundamental units of length and time are 
defined by /ho = yhjmujz and tho = 1/^^ respectively, 
with the dimensionless interaction strength being = 

m/TC'ujz/ (27r/^), Ir is the size of the BEC in the radial 
direction. 

To evaluate number of atoms in wells two and three in 
ID we use the definitions. 



N2{t)= / \^{z,t)\^dz (11) 
N3{t)= / \^{z,t)fdz. (12) 

Jzo 

In 3D, Eqs. (11) and (12) are generalized to 
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which are used to define the fidelities of the process shown 
in Figs. 1(b) and 4(a). 
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